Magnetic field driven instability of charged center in graphene 
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It is shown that a magnetic field dramatically afltects the problem of supercritical charge in 
graphene making any charge in gapless theory supercritical. The cases of radially symmetric poten- 
tial well and Coulomb center in an homogeneous magnetic field are considered. The local density 
of states and polarization charge density are calculated in the first order of perturbation theory. 
It is argued that the magnetically induced instability of the supercritical Coulomb center can be 
considered as a quantum mechanical counterpart of the magnetic catalysis phenomenon in graphene. 
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I. INTRODUCTION 

Recently it was shown that atomic collapse in a 
strong Coulomb field [H, H | , a fundamental quantum rel- 
ativistic phenomenon still inaccessible in high-energy ex- 
periments, can be readily investigated in graphene. In 
quantum electrodynamics (QED), taking into account 
the finite size of nucleus Q, theoretical works on the 
Dirac-Kepler problem showed that for atoms with nu- 
clear charge in excess of Z > 170 the electron states dive 
into the lower continuum leading to positron emission 
ii- 

In graphene, the effective Coulomb coupling constant 



is given by /? = Za/n, where a 



'-/hvF 



2.19 is 

the "fine-structure" coupling constant, vp ~ lO^m/s 
is the velocity of Dirac quasiparticles, and k is a di- 
electric constant. The Hamiltonian of the system is 
not self-adjoint when the coupling (3 exceeds the crit- 
ical value Pc = 1/2 [1-4]. Similar to the Dirac equa- 
tion in QED, one should replace the singular 1/r po- 
tential by a regularized potential which takes into ac- 
count the finite size of the charged impurity, R: V{r) = 

— ^0{t — R) — ^9{R — r). For gapped quasiparticles 
in such a regularized potential, the critical coupling is 
determined by /3c = 1/2 + ■k'^ / \og^{c/:^R/hvF) [1, where 
A is a quasiparticlc gap and the constant c 0.21, and 
Pc tends to 1/2 for A ^ or i? ~> 0. 

Since the electrons and holes strongly interact b" 
means of the Coulomb interaction, one may expect ^ 
[To| an excitonic instability in graphene with subsequent 
phase transition to a phase with gapped quasiparti- 
cles that may turn graphene into an insulator. This 
semimetal-insulator transition in graphene is actively 
studied in the literature, where numerical simulations 
give the critical coupling constant ac « 1-19 flll.[l^. 

In a many body system or quantum field theory, the su- 
percritical coupling leads to more dramatic consequences 
compared to the case of the Dirac equation for the elec- 
tron in the Coulomb potential. Unlike the case of the 
Coulomb center, the many body supercritical coupling 
instability cannot be resolved through a spontaneous cre- 
ation of a finite number of electron-positron pairs. Like 
the Cooper instability in the theory of superconductiv- 
ity, the QED supercritical coupling instability is resolved 
only through the formation of a condensate of electron- 
positron pairs generating a mass gap in the spectrum 



The presence of a magnetic field makes the situation 
even more interesting. It was shown in [l^ that magnetic 
field catalyses the gap generation for gapless fermions in 
relativistic-like systems and even the weakest attraction 
leads to the formation of a symmetry breaking conden- 
sate. Therefore, the system is always in the supercritical 
regime once there is an attractive interaction. The mag- 
netic catalysis plays an irnpqrtant role in quantum Hall 
effect studies in graphene [i3-[2iJl , where it is responsible 
for lifting the degeneracy of the Landau levels. 

The magnetic catalysis phenomenon suggests that the 
Coulomb impurity in a magnetic field in graphene should 
be supercritical for any Z . The Dirac equation for quasi- 
particles in graphene in the Coulomb potential in a mag- 
netic field was considered in [2l| where exact solutions 
were found for certain values of magnetic field, however, 
no instability or resonance was found. 

In QED in (3+1) dimensions, the Coulomb center 
problem in a magnetic field was studied in . There it 
was found that magnetic field B confines the transverse 
electronic motion and the electron in the magnetic field is 
closer to the nucleus than in the free atom. Thus, it feels 
stronger Coulomb field. Therefore, Zca decreases with 
B. This result is consistent with the magnetic cataly- 
sis phenomenon |14| . according to which, magnetic field 
catalyses gap generation and leads to zero critical cou- 
pling constant in both (3+1)- and (2+l)-dimensional the- 
ories. 

We would like to stress that the presence of an homoge- 
neous magnetic field changes qualitatively the supercrit- 
ical Coulomb center problem. Indeed, if magnetic field 
is absent, then the supercritical Coulomb center insta- 
bility leads to a resonance which describes an outgoing 
positron propagating freely to infinity. However, since 
charged particles confined to a plane do not propagate 
freely to infinity in a magnetic field, such a behavior is 
impossible for the in-plane Coulomb center problem in a 
magnetic field. Therefore, a priori it is not clear how the 
instability suggested by the magnetic catalysis manifests 
itself in the Coulomb center problem in a magnetic field. 
To answer this question is the main aim of this paper. 

In Sees, ini mil we consider the Dirac equation for the 
electron in the potential well and Coulomb center in a 
magnetic field. We study the local density of states 
(LDOS) and induced charge density for both cases in 
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Sec. IIVI where similarities and differences between the 
cases of gapped and gapless quasiparticles as well as the 
potential well and Coulomb interactions are discussed. In 
Sec. |V] we give a brief summary of our results. Finally, 
we provide the details of our calculations of the LDOS 
and polarization charge density in Appendix A. 



with j = ±1/2, ±3/2, . . . . For functions /(r), g{r) we 
get the following equations: 



II. POTENTIAL WELL 



, , J - 1/2 ^ r E-^A- V{r) ^ 
5+^— .9 + ^.9 J- / = 0, (2.6) 



The electron quasiparticle states in vicinity of the K± 
points of graphene in the field of Coulomb impurity and 
in an homogeneous magnetic field perpendicular to the 
plane of graphene are described by the Dirac Hamiltonian 
in 2+1 dimensions 



H = hvpTp + ^Ara + V{r) 



(2.1) 



where the canonical momentum p = — iV + eA/c in- 
cludes the vector potential A corresponding to the exter- 
nal magnetic field B, Ti are the Pauli matrices, and A is a 
quasiparticle gap. The two component spinor 5*^3 carries 
the valley (^ = ±) and spin (s = ±) indices. We will use 
the standard convention: '^^^ = {'4'a,iPb)k+s whereas 
^-s = (V'-B) ''I^a)k-s, and A, B refer to two sublattices of 
hexagonal graphene lattice. Since the interaction V{r) 
does not depend on spin, in what follows we will omit 
the spin index s. 

It is instructive to consider first the Dirac equation for 
the electron in a potential well V{r) = —VoO{rQ — r) with 
Vb > in a magnetic field perpendicular to the plane of 
graphene. We have 



hvp{- 



{E - V{r)) *(r) 



(2.2) 



where Di = di + {ie/hc)Ai with i = x,y is the co- 
variant derivative and the symmetric gauge {Ax, Ay) = 
{B/2){—y,x) is used for the magnetic field. It is clear 
that the solution at point is obtained from the solu- 
tion at point changing A — ^ — A and exchanging the 
spinor components ipA iJjb- 
In polar coordinates 



iDx + Dy 



iDx - Dy 



.d Id 
dr r 86 



ld_ 

r dcf) 



.d_ 

dr 



ieBr\ 
ieBr^ 



2hc 



(2.3) 



where I = y/hc/\eB\ is the magnetic length. These equa- 
tions are easily solved for the potential well in two regions 
r < Tq and r > Tq in terms of confluent hypergcomctric 
functions. In the region r < rg, eliminating the function 
g{r), we obtain the second order differential equation for 
the function /(r): 



f"--f' + 
P 



2pt 



1 _ f - 3 - 3/4 _ 

2 p2 4 



/ = 0, 



(2.7) 

and in the region r > Tq we have the same equation but 
with Vo = 0. Here we introduced the following dimcn- 
sionless quantities: 



2 _ P[{E + Vor-A^] , 
Pv ~ ;t77 — V? ' P 



2{hvF) 



P{E^~A^ 



(2.8) 



and p ~ r/l. We get the following solutions which are 
regular at ?- 0, 



hip) 



r(j + 1/2) 



T., ■ 1 2 . 1 



(2.9) 



9i{p) = 



l{E + Vo-(A) 



V2h- 



pi 



V{j + 3/2) 



^ , ■ 1 2 ■ 3 p2 



(2.10) 



and decrease at infinity, 

/2(p) = C2pi+"-e- 



pV4 



T I ■ 1 2 ■ 1 



52 (p) 



2 0+3 -0^/4 



V2hvFC: 
1{E + ^A) 

T r 1 2 . 3 p2 



(2.11) 



(2.12) 



We can represent 4'(r) in terms of the eigenfunctions of 
the conserved angular momentum = + (yyj2 = 
—id/d4> + cFz/2 as follows: 



respectively [note that these expressions are valid at all 
J =±1/2, ±3/2,...]. 

Then sewing solutions at r = tq. 



■h)4>f{r) 
9{r) 



(2.4) 



h{p) 



hip) 



giif) 



P=PO 



52 (p) 



P=Po 



Po = -r 



(2.13) 
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we obtain the following transcendental equation for en- 
ergies of solutions with the total angular momentum j: 



P{E + Vo- eA)a> J + i - pI,J + |; f 



(2.14) 



Below we analyze this equation analytically and numer- 
ically. 



A. Instability in the absence of magnetic field 

In this subsection, we consider the problem of the po- 
tential well instability in the absence of magnetic field 
{B = 0) which will serve us as a useful reference point 
in the next section where we study instability in a mag- 
netic field. For B = 0, the energy spectrum may be 
obtained either solving the Dirac equation from the very 
beginning, or taking the limit of zero field (l — > oo) in 
Eq. (|2.14p . In the last case one needs to use the following 
formulas [23| . 

$(a,c; z) = e''$(c - a, c; -z), (2.15) 
lim $(a,c; -z/a) = r(c)z^ Jc_l(2^/z), (2.16) 

a— >-oo 

lim [r(l + a - c)^'(a, c; -z/a)] 

a^oo 



i^,^i\(2^/i), lmz>0, (2.17) 



and |arga| < tt for the last two equations. 
Assuming \E\ < A we obtain 



E + Vo-^A J 



J + l/2 



Ve^ 



where J,(z), hI^\z) are the Bessel and Hankel 
functions, respectively, f3 = y/ {E + Vb)^ — A'^/hvp, 
(3' = V-B^ — A'^ /hvF, and the square roots are de- 
fined as Im/3,Im/3' > 0. In the regions with 
Im/3,Im/3' ^ one can use the relations H'i}\iz) ~ 
{2/TTi)e-^'"'/^K^{z),J^{iz) = e'^'^/^I^iz). Eq.d^TH]) is 
invariant under the change j — > — j, ^ — > — ^ |24| . 

Taking for the definiteness the point = — ), one 
can see from Eg. p.lSp that the energy spectrum is con- 
tinuous for \E\ > A and a discrete one for \E\ < A. The 
first bound state E < A appears at an arbitrary small in- 
teraction Vq. Indeed, taking j — —1/2 that corresponds 
to the smallest centrifugal barrier, we find 



(1) 

J + l/2 



(/?Vo) ' 



(2.18) 



~ A 



1 



-27 



(2.19) 



where 7 is the Euler constant. Note that there is no 
solution with the energy £^ < A at the K- point with 
angular momentum j = 1/2, but such a solution exists 
at the i^-L point similarly to the case of the Coulomb 
potential Q- 

As Vb grows, at the critical strength of interaction 



Vr 



cr 




(2.20) 



(jo,i ~ 2.41 is the first zero of the Bessel function Jq{x)) 
the lowest in energy bound state dives into the lower 
continuum {E ~ —A). We note that for the zero gap 
case (A = 0) there are no bound states at all. In the 
supercritical regime for Vb > Vbcr = f>'VFjo,i/''^o (which 
follows from Eq. (|2.20p at A = 0) resonances with com- 
plex energies appear leading to instability of the system 
similar to the case of the supercritical Coulomb center 
[l], @ . The occurrence of resonant states synchronously 
with diving into the lower continuum of the lowest in en- 
ergy bound state is the standard characteristic of QED 
systems [H, 0] . We will see in the next subsection that 
the presence of an homogeneous magnetic field changes 
this conclusion. 

Near the critical value of coupling the energy of reso- 
nant state is given by 



E 



5 = 



Vb — Vbcr 



■ exp 



ln(l/<5) V21n(l/<5)j' 
{Vo - l/ocr)roe^-i 



2Hvp 



0<(5<1. (2.21) 



The dependence of energy on the Vq — Vbcr (deviation 
from the critical value) for the potential well interaction 
is nonanalytical one but differs from the essential singu- 
larity that takes place in the Coulomb center problem. 
This, of course, is related to the absence of scale invari- 
ance for the potential well V{r). 



B. Magnetically driven instability 

Before we consider the instability of the potential well 
problem in an homogeneous magnetic field it is useful 
to recall the Landau energy levels for the electron states 
in graphene in a magnetic field. If the interaction van- 
ishes (Vb = 0,ro — >■ 0), Ea. (|2.14p gives the well known 
spectrum of Landau levels: 



E : 



<A, j< 



1 

'r 



(2.22) 



£; = ±WA2 + 2n 



fhvFV 1 

— j ,n = l,2,...,, + -<n. 



(2.23) 

Note that the level E ~ A {E = —A) is present only at 
the {K+) point. 
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For nonzero Vq , the Landau energy levels are no longer 
degenerate. Using the sewing equation (|2.14p . we can 
determine the evolution of degenerate solutions with Vq . 
For solutions of the Landau level E = A with different 
j, their energies as function of Vq/A (at fixed magnetic 
field B) are plotted in Fig. [U for lA/{y/2hvF) = 0.1 and 
Po = ''o/^ = 0.02. We sec that as Vq increases more 



-1/2 Eq. ([2?T4)) becomes 




FIG. 1: The evolution of degenerate solutions of the lowest 
Landau level at the K- point as a function of the dimension- 
less ratio Vo/A. 

and more solutions with different j cross the energy level 
E = — A. As usual [1, Q, this means that vacuum of the 
second quantized theory is unstable with respect to the 
creation of electron-hole pairs. However, as we discussed 
in Introduction, there are no resonance solutions in the 
presence of constant magnetic field. The reason for that 
is the presence of the positive r^/4P term due to the 
magnetic field in the effective Schrodinger-like equation 
for one component of the spinor function (see Ea. p.4p 
and Fig. 3 in the next section) which qualitatively changes 
the asymptotic of the effective potential: in nonzero field 
it grows at infinity instead of decreasing as in the case 
_B = 0. Therefore quasiparticles are confined in such 
a potential and cannot escape to infinity forming only 
discrete levels. 

We would like to note that the situation under consid- 
eration is analogous to that for a deep level vacancy in a 
many electron atom. There electron states as solutions 
of the Dirac equation in the Coulomb potential of the nu- 
cleus are stable. However, taking into account the inter- 
action with the second quantized electromagnetic field, 
the electrons on higher energy levels are unstable with 
respect to the transition to the vacant state with emis- 
sion of photons. 

The critical potential Vocr is defined as the potential for 
which the first crossing occurs. According to Fig.l, such 
a crossing is first realized for the state with j = —1/2 
(the potential well interaction lifts the degeneracy of the 
Landau levels in quantum number j). Let us analyze in 
detail how this state evolves with Vq. For the state with 



2 <I>(-p^,l;f 



{E + Vq- A) 



= -{E-A) 



where we used the relation 

lim ^^"'"^'^^ = r(a-H7H-l) 
c->-m r(c) r(a)(m + l)! 

x<f>(a + m + l,m-|-2;a;), m = 0,l, 





'-p',0; 


i) 






2 ) 



(2.24) 



(2.25) 



For Vq ^ 0, Eq. (|2.24[ ) implies the following bound state 
at the K_ point: 



E 



Voil 



(2.26) 



that is in contrast with nonanalytical behavior in the 
coupling Vq in the absence of magnetic field described by 
Eq. (|2.19p . [At the A'+ point a similar bound state exists 
but with angular momentum j = +1/2.] 

As the coupling Vq grows, energy of this bound state 
decreases and finally crosses the level E = —A at some 
critical value Vbcr- For E = —A, p'^ = 0, py = 1^{Vq - 
2AVo)l2{hvFf, using 'J'(0,0;z) = 1'(0,l;z) = 1 we find 
that Ea. (|2.24p defines the following equation for Vocr'- 



2A 



2$ 



a,l,f 



p2$ l + a,2, 



(2.27) 



where a = -/ Vocr(Vbcr - 2A)/2{hvFY. [Note that 
at zero magnetic field (/ = oo) Eq. (|2.27p reduces to 
Eq. (|2.20p for Vbcr that tends to a finite value in the gap- 
less limit.] The critical potential strength Vbcr as a func- 
tion of A is plotted in Fig. [5] for different values of the 
parameter po which defines the ratio of the potential well 
width to the magnetic length. Analytically, it is not dif- 
ficuh to find that for po < 1 Eq. (P??f|) implies 



1/oc,. = 2A(1 + 2P/rl) 



(2.28) 



It is clearly seen from Figl^l and from Eq. (|2.28p that the 
critical potential strength Vbcr decreases with the growth 
of a magnetic field (or, with the decrease of /) at fixed j-q 
and A. The physical reason for that is that the magnetic 
field forces electron orbits to become closer to the charge 
center, making attraction stronger and, thus, effectively 
lowering the critical coupling. 

What is surprising here is that Vbcr tends to zero as 
A — > 0. Thus, the presence of an homogeneous magnetic 
field leads to the instability of gapless quasiparticles in 
the second quantized theory for any value of the poten- 
tial strength Vb- This result suggests that the Coulomb 
center in gapless graphcne in a magnetic field may be also 
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rp Vqc 



means of the relation 





Po=10 




Po = l 




p„=0.5 




p„=0.2 




p„=0 


^ 



0.0 



0.5 



1.0 
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2.0 



FIG. 2: The critical potential Vbcr as a function of a gap 
for different values of po- The case of zero magnetic field 
corresponds to po = 0. 



unstable for any value Ze, the problem which we study 
in the next section. 

Finally, we will analyze states with energies near ±A 
and large by modulus negative momenta j. We find that 
there exists an infinite series of levels approaching the 
energies ±A asymptotically at large \j + 1/2| (i.e., for 
sufRciently large j the effect of the potential interaction 
Vo can be neglected and the Landau levels are recovered) . 
For Vq — >■ 0, they behave as 



k+1 



» 1. 



r(A: + 1) V 2 

(2.29) 

This can be found directly by solving Eq. (|2.14p . first tak- 
ing there the limit j + 1/2^-— /cby means of Eq. ()2.25|) 
and analyzing then the equation at weak coupling and 
large k. Alternatively, Eq. (|2.29P is obtained as the first 
order correction in the interaction to the levels ±A at K± 
points in a magnetic field. Note that the levels (|2.29p lie 
below A for the point and below —A for the Kj^ 
point, respectively. 



we get the Schrddinger-like equation, 

- X"ir) + U{r)x{r) = £x{r) 

where 

E = E^ - A^ 
and the effective potential, U = Ui 
Vi2E^V) , jU + 1) 



(3.1) 



(3.2) 



(3.3) 



Ui 



4^4 



J - 1/2 



,(3.4) 



V" 



r 

2P 



V 2 
2V' 



V 



E-^A-V 



E-^A-V 



(3.5) 



We plot effective potential U{r) near point for E = 
—A and j — —1/2 in Fig. 3. There the energy barrier 
in the absence of magnetic field is clearly seen, which 
leads to the appearance of resonances for sufficiently large 
charge. The presence of non-zero magnetic field changes 
the asymptotic of the effective potential at infinity and, 
thus, forbids the occurrence of resonance states. This fea- 
ture distinguishes qualitatively the Coulomb center prob- 
lem (as well as the potential well problem) in a magnetic 
form from that at i3 = 0. 



U(r) 




FIG. 3: The potential U(r) as a function of a distance from 
the Coulomb center at zero and nonzero magnetic field for the 
state with E = —A and j — —1/2. 



III. THE COULOMB CENTER 

The equations for the functions /(r) and g{r) for the 
Coulomb center problem follow directly from Eas. (|2.5p 
and (Hill) by setting V{r) = -Ze^/r e{r - i?) - 
Ze^ jR 9{R — r) there (we take the dielectric constant 
K ~ 1). Eliminating, for example the function /(r). 
one can get a second order differential equation for the 
function g{r). Further, introducing the function x('') by 



Unfortunately, Ea. p.2p belongs to the class of equa- 
tions with two regular and one irregular singular (at 
r ~ oo) points, and exact solutions of this equation can- 
not be expressed in closed form in terms of known special 
functions. 

Since we are interested in solving Eas. (|2.5p , (|2.6p with 
the Coulomb potential in the regime Za — > (a = 
e'^/hvp), we can find solutions using perturbation the- 
ory. For Za = 0, solutions of Eqs. (j2.5p . (|2.6p are the 
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well known that Landau states degenerate in the total 
angular momentum j. For the level E^^') = A their nor- 
malized wave functions ()2.4p have the form (at the A'_ 
point) 




(3.6) 



where k = -{j + 1/2) = 0,1,2,.... The Coulomb po- 
tential removes degeneracy in j. Energy corrections of 
perturbed states of the Landau level E'-^^ = A are found 
from the secular equation 

Since V^jfe^ is a diagonal matrix, we easily obtain 



El 



Vkk = - 



k\2H 



dp p''" e 



2k „-pV2 



Ze^r(fc+1) 



(3.7) 



lV2T{k + 1) ■ 

Thus at large k the energy levels accumulate near E = A: 



Ek-A 



W2k 



(3.8) 



Like in the case of the potential well interaction consid- 
ered in the previous section, the largest correction by 
modulus Sq^' = —Zahvpy/Tr/lV^ occurs for the state 
with j ~ —1/2 {k = 0). The critical charge is deter- 
mined by the condition E ~ E^^") + E^'^^ = — A when the 
level E crosses the level of filled state. This gives 



Zra 



2V2AI 



(3.9) 



Like in the case of the potential well in a magnetic field, 
the critical charge p.9p tends to zero as A 0. This 
means that magnetic field indeed dramatically affects the 
Coulomb center problem in graphene making any charge 
in gapless theory supercritical. Eg. p.Op gives the crit- 
ical Coulomb coupling in the regime Za — > in the 
first order of perturbation theory. For arbitrary values 
of Za, we calculated the dependence of the critical cou- 
pling on the gap numerically. The corresponding results 
arc presented in Fig. |4l where, for the parameter regu- 
larizing the Coulomb potential, we took R = 10~'^L The 
dashed red line in Fig. 2] gives the critical Coulomb cou- 
pling ZcCt = 1/2 + -K^ / \o^{cAR/hvF) in the absence of 
magnetic field (see Fig. 1 in Ref.Q). Thus, at weak 
magnetic field [l — > 00) the critical coupling tends to 
1/2 while ZcCt for ZA in the gapless or strong 
magnetic field regime. 




0.0 



0.2 



0.4 
lA/^iVf 



0.6 



0.8 



FIG. 4: The critical Coulomb coupling ZcCi as a function 
of gap at zero (dashed red line) and nonzero magnetic field 
(dotted black line) for the state with j = —1/2. The straight 
black line corresponds to the critical Coulomb coupling in the 
first order of perturbation theory given by Ea. (|3.9p . 



IV. THE LOCAL DENSITY OF STATES 

It is interesting to see how a magnetic field and the 
Coulomb center affect the local density of states of quasi- 
particles that can be directly measured in scanning tun- 
neling microscope (STM) experiments. The crucial dif- 
ference of the case of gapless quasiparticles from that of 
gapped ones in a magnetic field is that the critical charge 
is zero for gapless quasiparticles, therefore, energies of all 
previously degenerate states of the lowest Landau level 
become negative. 

The LDOS at the distance r from impurity is given by 

p{E;y) = --iTlmG{Y,Y;E + ivi), t] ^ 0, (4.1) 

where trace includes the summation over valley, sublat- 
tice and spin degrees of freedom. The retarded Green's 
function G(r, r';E + irj) in a constant magnetic field can 
be written in the form 

G(r, r';E) = e**('-'"-')G(r, r'; E), (4.2) 

r 

'fir,r') = Ar\z)dz\ (4.3) 

r' 

where we separated the gauge dependent (Schwingcr) 
phase $(r, r') from a gauge invariant part of the Green's 
function G{r,r';E) . The last one satisfies the following 
Lippmann-Schwinger equation: 

G(r, r';E) = Go{r-r';E) + J dr"Go(r - r"; E) 

xy(r")G(r",r';i;)e'l*('''"'")+*^''"'""')+*('"'^'")l.(4.4) 

[Note that the Green function G{r,r';E) is not transla- 
tion invariant in presence of an impurity unlike the non- 
interacting function Go(r — r'; E).] For weak interaction 
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we can calculate the LDOS in the first order in the per- 
turbation theory, 



p{E-v)= p^{E-v)+5p{E;v), 



(4.5) 



where po(£';r) is the LDOS for free quasiparticles in a 
magnetic field, and 

5p{E- r) = -;^Im j dr'tr Go(r - y')V{y')Gq{y' - r) 

(4.6) 

First we consider the case of gapless quasiparticles. In 
this case, po{E] r) and 5p{E] r) are calculated for gapless 
quasiparticles in Appendix. We got there that the LDOS 
decreases at large distances (r >> ro,/) as 



5p{v,E) 



2T:{hvp)' 



:Im[Ar2(-A)]e 



2/2 



2A 



(4.7) 



for the potential well, while for the Coulomb center we 
obtained 



Spir,r,E) 



1 



K 2T:{hvp) 



2r3 



(4.8) 



with the functions Bi{X) given by Eqs. (|6.14p . (|6.15p and A 
is defined after Ea. (|6.4p . For the induced charge density. 



dE5p{v,E), 



(4.9) 



using Eas. (|4.7p and (j4.8p . we find that it is positive at 
large distances and decreases exponentially for the po- 
tential well and, due to j'^^dEBo{\) = 0, as l/r^ for 
the Coulomb interaction. 



ZeH 1 
nhvp 



3C(-l/2) 
27rV^ 



0.07. (4.10) 



We remind that in the absence of magnetic field the po- 
larization charge diminishes as l/r^ both in the super- 
critical Zajn > 1/2 Coulomb center [l| and potential 
well Vb > Vbcr = fiVFjo.i/fo interactions. 

The situation is quite different in the case of gapped 
quasiparticles. Here we will consider the polarization 
charge density for the most interesting case of the 
Coulomb center in an homogeneous magnetic field. The 
polarization charge density (|4.9p can be rewritten in the 
more familiar form 



riindir) 



_E<-A 



^^(r) ^\^JJE{r)\' 



(4.11) 



where ipp and are the Landau wave functions and 
the wave functions of the Coulomb center problem in a 
magnetic field, respectively. Since we consider the case 
where Za is small, the corrections to the wave functions 
of negative energy states of the deep Landau levels de- 
fined by Eq. (|2.23p can be ignored. We will consider only 



the corrections to the lowest Landau level states given by 
Eq. (11221) • 

In the first order of perturbation theory, wave functions 
of the Landau level E = —A are sought as superposi- 
tion of all degenerate states with unknown coefficients 

E = Sj<i/2 '^i'^'^-A' 'wliGre ip-A wave functions of 
the Landau level with E = — A and total momentum j. 
The unknown Cj are determined by the equation 



3132/ '-32 



0. 



(4.12) 



32<-l/2 



Since 



3132 



^3131^3132 is a diagonal matrix, the sec- 



ular equation \Vj^j^ — iS^^^j = is trivially solved and 

E'j^' = Vjj. Then Eq. (|4.12p implies that the correspond- 
ing unknown coefficients Cj-^ are equal Cj^ = 1 for ji = j 
and zero for all other ji . Notice that Cj = 1 because the 
wave functions arc normalized. This means that 

the wave functions of perturbed states of the Landau 
level £^ = — A do not change in the first order of per- 
turbation theory. Consequently, according to Ea. (|4.1ip . 
they do not contribute to the polarization charge den- 
sity. Clearly, the polarization charge appears only when 
the first perturbed state of the Landau level E = A with 
j = —1/2 crosses the threshold of filled states of the low- 
est Landau level £: = -A. Using Ea. (|X^ . (gUI]), and 
the fact that perturbed wave functions of the E ~ A 
Landau level states do not change in the first order in 
the Coulomb potential, similarly to the case of perturbed 
wave functions of the E = —A Landau level considered 
above, we conclude that for the critical charge ZcO: given 
by Eq. (|3.9p the polarization charge density equals 



ni„d{r) 



2ttP 



(4.13) 



Thus, the polarization charge density is concentrated 
near the impurity where it is negative and quickly de- 
creases at large distances. 



V. CONCLUSION 

In this paper we showed that in an external magnetic 
field the value of the critical coupling for the onset of 
instability of a system of planar Dirac gapless quasipar- 
ticles interacting with charged center (charged impurity) 
reduces to zero. This result serves as a quantum mechan- 
ical single-particle counterpart of the magnetic catalysis 
phenomenon in graphenc. The cases of radially sym- 
metric potential well and Coulomb center were analyti- 
cally considered. The local density of states and induced 
charge density were calculated in the first order of per- 
turbation theory for gapless quasiparticles. 

The crucial ingredient for the instability is the exis- 
tence of zero energy level for gapless Dirac fermions in a 
magnetic field which is infinitely degenerate. In this case 
any weak attraction leads to the appearance of empty 
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states in the Dirac sea of negative energy states and to 
the instabihty of a system. 

One should stress a quahtativc difference in the phe- 
nomenon of instabihty between gapped and gapless 
quasiparticles. In the case of gapped quasiparticles there 
is a finite critical value for the strength of interaction 
when the lowest unfilled level crosses the first filled one 
forming a hole in the sea of filled states. As the coupling 
grows, more and more levels cross that level. Clearly, 
the system tries to rearrange itself filling in empty states 
whose presence is a signal of instability. The important 
difference of the case of gapless quasiparticles from that 
of gapped ones, besides the critical coupling being zero, 
is that an infinite number of states of the previously de- 
generate lowest Landau level become vacant. 

Thus, the presence of an external magnetic field 
changes dramatically the problem of atomic collapse in 
graphene in a strong Coulomb field [l[ . Clearly the prob- 
lem becomes a many body one and requires field theo- 
retical methods to find a true ground state. One should 
expect that the gap generation for initially gapless quasi- 
particles should take place already in the weak coupling 
regime in the presence of a magnetic field [l6j . 
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can be explicitly performed by means of the formula 



E 

ri=0 



n + b 



= F(6)*(fe;l + a;x) 



(6.2) 



(see, Eq.(6.12.3) in the book [2J|), leading to a closed ex- 
pression for the free Green's function (see, recent papers 

MM), 



Go{r;E) 



4:TTh?Vp 



E 



P_F(-A)* -A;l; 



r 

2/2 



P+r(f-A)*(l-A;l;^ 



+ z?i,;;,^F(l_A)* (^I-A;2;i^j|. (6.3) 

Here \E'(a;c;x) is the confluent hypergeometric function 
which is related to the Whittaker function. 



c- 1 



^'(a;c;a;) = e^x ^ ^^W^,.^|,{x), n = - - a, ^. = 

(6.4) 

r(x) is the Euler gamma function, and A = (£' + 
irjfP/{2h^vl). 

The LDOS of free quasiparticles in a magnetic field 
does not depend on r and is given by 

po{E) = lim LTitr[Go(r; E + ii])] = J „ 



X hm Im < [E + iij) 



r(-A)* -A;l 



r 

2P 



-F(I-A)1' 1-A;1 



2P 



(6.5) 



VI. APPENDIX 



The Green's function of free quasiparticles in a mag- 
netic field is well known (see, for example, [13, and 
in the configuration space it has the form of a series over 
the Landau levels (we consider the zero gap case). 



Go{r;E) 



27r/2 



M2 









E 


(^) 


+ ^+^.-1 (^) 









(6.1) 



where M„ = {Uvp /l)y'2n are the energies of the Landau 
levels, F± = (1 ± T3)/2 being the projectors, L'^{z) the 
generalized Laguerre polynomials (by definition Lniz) = 
L'^{z) and L"^{z) = 0), and the Pauli matrices ti_2 act 
in the sublattice space. The sum over the Landau levels 



The hypergeometric function \['(a;c;a;) at small x be- 
haves as 

*(a;l;a:) ~ --^[Inx -I- V(a) + 27] + 0(a;lna;), 
r(a) 

*(«; 2; x) ~ — i— + — ^— [In x + V^a) 
i [ajx 1 (a — 1) 

+27- 1] + 0(xlna;),(6.6) 
where ipiz) is the digamma function. Therefore 
1 



Po{E) 



^lui[{E + i5) (^(-A)+^'(l- A))], 

(6.7) 

and the LDOS free quasiparticles in a magnetic field fi- 
nally is found to be 



P.[E) . -I 



S{E) + ^ [S{E - M„) + 5{E + M„)] 



(6.8) 

(compare with Eq.(4.2) in Ref.j27l|). 

The first order correction to the LDOS due to the in- 
teraction is given by Ea. (|4.6p . For the radial well to find 
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the asymptotic at distances r ^ rp , where Tq is a range of 
the potential, we can put r' = in the arguments of the 
free Green's functions in Eq. ()4.6p and get the following 
behavior: 



Spir,r;E) - Vor^olmtT[Goir; E)Go(-r; E)] 
2Vorl 



2{T:hvpl) 



2lm[AV'(-A)]ln— , 



,Im[Ar2(-A)]e-V2q^) , (6.10) 



,2 \ 2|A| 
2/2 



in the regions I >> r >> 7'o and r » max(/,ro), re- 
spectively. As is seen, the correction to the free LDOS is 
an odd function of energy. 

To find the first order correction due to the Coulomb 
potential we first calculate the correction to the Green's 
function which is given by 

<5G(r, r; E) = -— [ dr'Go(r - r'; E)j^Go{r' - r; E) 
K J \r'\ 

' dv'Go{v'-.E)y—^G^{-v'-E). (6.11) 

Taking trace over spin and Dirac indices, performing in- 
tegration over the angle by means of the formula 



27r 



die 



K 



Arr' 



v/r2 +r'2 - 2rr'cos6' r + r' " \ [r + r'Y j ' 

(6.12) 

where K{k) is the elliptic integral of the first kind 
{K{Q) = 7r/2), and calculating the imaginary part, we 
obtain 



5p{v,v,E) 



2Ze^ sgn{E) f dr'r' _^ 



K (jrhvpl)^ J r + r' 




K 



Arr' 



(r + r'Y 



^ [A([L„(a;)]2 + [L„_i(a;)] 



4-2x[i^_i(x)]2]<5'(A-n) 

OO 

- ^ [A(L„(a;)i,„(a;) +i„_i(a;)i,„_i(a;)) 

+2xLL,(.)X^_,(.)] ^(A-»)-^(A-^) l 

n — m J 

where x = r'^/2P. The correction to the LDOS at large 
distances r >> / is given by Ea. (|4.8p where the energy 



dependence is given by the functions 
Bo(A) =sgn(i;) 



(6.9) Bi{X) = sgn{E) 



A(S'(A)-l-2^(A + n),5'(A-7 

n=l 

CXD 



(6.14) 



Tl = l 



25{X - 1) + 2 ^[A(2?7, + 1) + 2n{n + 1) 



ri=l 



X ((5(A-?i-l)-(5(A-n))] 



(6.15) 



To calculate the integrals of the Laguerre polynomials we 
used the following generating function (see. Appendix A 
in [HI): 



/"™(y)= / dte-H"M2V^t)LZ{t)LfJt) 



ml 



e-yL^-^{y)L-^;i{y). (6.16) 



Expanding the left and right sides in y we find the stan- 
dard orthogonality relation. 



'dte-*rLg(t)L^(0^ (6.17) 



and the integrals 



dte-HL^{t)Lra{t) = (-1)"+™ [(2n+l),S„„ 
+ (m + l)(5,i,,„+i + (n + l)(5„i,„+i] , (6.18) 
dt e-H^L\{t)Ll{t) = (-l)"+"(m + l)(n + 1) 



(6.19) 



[2(5„m + 5n,m+l + 5m,n+l\ 



For the integrals of the functions Bi{X) we get 
/°^dSSo(A) =0, while 



j dSSi(A) = ^6V2C(-l/2), (6.20) 

— OO 

where C,{z) is the Riemann zeta function. 
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